The Astrophysical Journal, submitted 

Preprint typeset using WT^^ style cmulatcapj v. 6/22/04 



EFFECTS OF COOLING AND STAR FORMATION 
ON THE BARYON FRACTIONS IN CLUSTERS 

Andrey V. Kravtsov^'^'^, Daisuke Nagai^"', Alexey A. Vikhlinin'' 

The Astrophysical Journal, submitted 



in 
o 
o 

(N 



(N 

> 

(N 
(N 

O 

in 
o 

o 



ABSTRACT 

We study the effects of radiative cooling and galaxy formation on the baryon fractions in clusters 
using high-resolution cosmological simulations that resolve formation of cluster galaxies. The simula- 
tions of nine individual clusters spanning a decade in mass are performed with the shock-capturing 
eulerian adaptive mesh refinement 7V-body-|-gasdynamics ART code. For each cluster the simulations 
were done in the adiabatic regime (without dissipation) and with radiative cooling and several phys- 
ical processes critical to various aspects of galaxy formation: star formation, metal enrichment and 
stellar feedback. We show that radiative cooling of gas and associated star formation increase the total 
baryon fractions within radii as large as the virial radius. The effect is strongest within cluster cores, 
where the simulations with cooling have baryon fractions larger than the universal value, in contrast 
to the adiabatic simulations in which the fraction of baryons is substantially smaller than universal. 
At larger radii (r > J'soo) the cumulative baryon fractions in simulations with cooling are close to 
the universal value. The gas fractions in simulations with dissipation are reduced by « 20 — 40% 
at r < O.Srvir and w 10% at larger radii compared to the adiabatic runs, because a fraction of gas 
is converted into stars. There is an indication that gas fractions within different radii increase with 
increasing cluster mass as /gas oc M^^^. We find that the total baryon fraction within the cluster virial 
radius does not evolve with time in both adiabatic simulations and in simulations with cooling. The 
gas fractions in the latter decrease slightly from z = 1 to 2; = due to ongoing star formation. 

Finally, to evaluate systematic uncertainties in the baryon fraction in cosmological simulations we 
present a comparison of gas fractions in our adiabatic simulations to re-simulations of the same objects 
with the entropy-conserving SPH code Gadget. The cumulative gas fraction profiles in the two sets 
of simulations on average agree to better than « 3% outside the cluster core {r/r^-^ > 0.2), but differ 
by up to 10% at small radii. The differences are smaller than those found in previous comparisons of 
eulerian and SPH simulations. Nevertheless, they are systematic and have to be kept in mind when 
using gas fractions from cosmological simulations. 

Subject headings: cosmology: theory-clusters: formation- methods: numerical 



1. INTRODUCTION 

In hierarchical models of structure formation clusters 
of galaxies form via collapse of a large representative vol- 
ume and are therefore expected to contain the baryons 
and dark matter in proportions close to the universal 
average. This fact together with independent measure- 
ments of the universal baryon fraction can be used as 
a powerful constraint on the mean density of matter in 
the universe (e.g.. White et al. 1993 : Evrard 1997). The 
apparent evolution of baryon fraction with redshift can 
be used to put constraints on the energy density content 
of the universe and properties of dark energy fSASakj 
119961 |Penlll997j) . The universality of cluster baryon frac- 
tions can be used to estimate total cluster mass and 
construct mass functio ns of clusters at different redshifts 
()Vikhlinin et alJl200^ . 

Recently, deep X-ray observations of nearby and 
high-redshift clusters have been used to measure the 
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cluster baryon fraction and co nstrain the power spec- 
trum shape and normalization ijVoevodkin fc VikhlininI 
lin^), the matter and dark eiiergy density of the 
universe llAllen et alJl200l l2004t [Vikhlinin ct al. 200^ 
iChen &: Ra tra"2004"). and even the neutrino mass and 
the equation of state of dark energy L Allen et al.. 2004'). 
Sunyaev-Zeldovich (SZ) effect observations combined 
with X-ray observations have also been used to measure 
the cluster baryon fraction a nd constrain the m atter en- 
ergy density of the universe lIGrego et a,lJl2001l^ . These 
measurements are complementary to (and competitive 
with) cosmological constraints using supernovae type la 
and anisotropics of the cosmic microwave background 
and are likely to be improved. However, they rely on 
the key assumptions that the baryon fraction within a 
given radius can be measured reliably, is independent of 
redshift, and can be converted to the universal fraction 
by the correction factor derived from cosmological simu- 
lations. 

In order to improve the baryon fraction based cos- 
mological constraints and to better understand associ- 
ated systematic uncertainties, these assumptions have 
to be carefully tested with simulations. During the 
last decade, a number of studies have considered the 
cluster baryon fraction and its evolution in cosmolog- 
ical simulations in the adiabatic regime (i.e., without 
radiative cooling). The main results of these stud- 
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ies are: (1) the baryon fraction within the cluster 
virial radius evolves only very weakly with time and 

(2) is in gener al » 10 — 1 5 % lower than the uni- 
versal average llEvrardI llQQOt iMetzler &: EvrardI 19941: 
Navarro et"^| 1995': Lubin et al." 1996': 'Eke et al."199??: 
Frenk et al.lIT999 : Mohr et al 1999: Bialck et al...2nOU : 

(3) this "baryon dep letion" can b e enhance d by strong 
pre-heating of gas fBiale k et al.l 1200 li iBorgani et al.l 
12002: Muanwong ct al. 2002t IKav eTall I2003D Note, 
however, that high-resolution eulerian simulations of 
lAnninos fc Norman ( 1996 ) exhibit baryon fraction larger 
than universal. Moreover, systematic comparison of adi- 
abatic simulations of a r nassive cluster pe rformed with 
different numerical codes ijFrenk et al.ll999l) has revealed 
systematic differences. The eulerian shock-capturing 
codes have generally produced clusters with little or no 
"baryon depletion," while gas fractions in SPH simula- 
tions were « 10% below the universal value. This in- 
dicates that there are systematic differences in gas and 
dark matter density profiles in simulations done with dif- 
ferent numerical schemes. 

More recently, several groups have presented analysis 
of baryon fractions in simulations that in clude ga s cool- 
ing and star formation llMuanwong et al.l2 002: Kav et al"! 
l2003t iVald arnini" 20031 lEttori et alJl200^ . These simu- 
lations have shown that cooling can increase the total 
baryon fraction (gas -I- stars) and change its dependence 
on cluster mass and cluster-centric radius compared to 
the adiabatic simulations. 

The purpose of the present paper is twofold. First, 
given the sy stematic differences between numerical codes 
identified in lFrenk et al.l l)1999(l . we would like to present 
a systematic study of the cluster baryon fractions in eu- 
lerian adaptive mesh refinement (AMR) simulations and 
compare them to the results obtained with modern SPH 
codes. Second, we investigate the effects of radiative 
cooling and star formation on the baryon fractions in 
the very high-resolution AMR simulations that resolve 
formation of individual cluster galaxies. 

The paper is organized as follows. In §[21 we describe 
our simulations and implementation of radiative cooling 
and star formation. In §|21we compare gas fraction pro- 
files in our adiabatic simulations with the simulations of 
the same clusters performed with the entropy-conserving 
Gadget code. In §^we discuss convergence tests and as- 
sess effects of resolution. In § [S] we present our main 
results and compare them to the previous results in § El 
Finally, in §[3we discuss our conclusions and implications 
of the presented results for cluster observations and their 
use as cosmological probes. 

2. SIMULATIONS 

In this study, we analyze high-resolution cosmologi- 
cal simulations of nine cluster-size systems in the "con- 
cordance" flat ACDM model: fl^ = 1 - ^Ia = 0.3, 
f]b = 0.04286, h = 0.7 and as = 0.9, where the Hub- 
ble constant is defined as 100ft. km s~^ Mpc~^, and as 
is the power spectrum normalization on 8h~^ Mpc scale. 
The simulations were done with the Adap tive Refinement 
Tree (ART) iV-body-|-gasdynamics code l)Kravtsovll999t 
^ravtsov et al. 2002) , a Eulerian code that uses adaptive 
refinement in space and time, and (non-adaptive) refine- 
ment in mass (Klvpin et alJl200l| to reach the high dy- 
namic range required to resolve cores of halos formed in 



self-consistent cosmological simulations. 

To set up initial conditions we first ran a low resolu- 
tion simulation of 8Qh~^ Mpc and 120ft.~^ Mpc boxes 
and selected nine clusters with the virial masses ranging 
from w 7 X lO^^ ^ 3 2 x lO^/i-^ Mq. The properties 
of clusters at the present epoch are listed in Table ^ 
The perturbation modes in the lagrangian region corre- 
sponding to the sphere of several virial radii around each 
cluster at z = was then re-sampled at the initial red- 
shift, Zi = 49 for eight clusters and Zi = 25 for the most 
massive cluster in the sample. For the Coma-size clus- 
ter we have resampled radius of l.5Rvir{z = 0), while 
for the rest of the clusters the resampling sphere had 
radius of 5i?vir- During the resampling we retained the 
previous large-scale waves intact but included additional 
small-scale waves, as described bv iKlvpin et al. (2001). 
The resampled lagrangian region of each cluster was then 
re-simulated with high dynamic range. 

High- resolution simulations were run using 128"^ uni- 
form grid and 8 levels of mesh refinement in the computa- 
tional boxes of 80/i^^ Mpc for CL2-CL9 and 120h^^ Mpc 
for the Coma-size CLl, which corresponds to the dy- 
namic range of 128 x 2* = 32768 and peak formal res- 
olution of 80/32,768 ~ 2AAh~^ kpc, corresponding to 
the actual resolution of « 2 x 2.44 « 5ft.~^ kpc. Only 
the region of ~ 3 — lO/i^^ Mpc around the cluster was 
adaptively refined, the rest of the volume was followed on 
the uniform 128'^ grid. The mass resolution corresponds 
to the effective 512^ particles in the entire box, or the 
Nyquist wavelength of Any — 0.469/i~^ and 0.312/i~^ 
comoving megaparsec for CLI and CL2-9, respectively, 
or 0.018/i~^ and 0.006/i^^ Mpc in the physical units at 
the initial redshift of the simulations. The dark mat- 
ter particle mass in the region around the cluster was 
2.7 X lO^ft-i Mq for CL2-CL9 and 9.1 x lO^/i"! Mq for 
CLl, while other regions were simulated with lower mass 
resolution. 

As the zeroth-level fixed grid consisted of only 128^ 
cells, we started the simulation already pre-refined to the 
2nd level {I = 0,1,2) in the high-resolution lagrangian 
regions of clusters. This is done to ensure that the cell 
size is equal to the mean interparticle separation and all 
fluctuations present in the initial conditions are evolved 
properly. During the simulation, the refinements were 
allowed to the maximum I = 8 level and refinement cri- 
teria were based on the local mass of DM and gas in each 
cell. The logic is to keep the mass per cell approximately 
constant so that the refinements are introduced to fol- 
low the collapse of matter in a quasi-lagrangian fashion. 
For the DM, we refine the cell if it contains more than 
two dark matter particles of the highest mass resolution 
specie. For gas, we allow the mesh refinement, if the cell 
contains gas mass larger than four times the DM particle 
mass scaled by the baryon fraction. In other words, the 
mesh is refined if the cell contains the DM mass larger 
than 2(1 — fb)fnp or the gas mass larger than = 4/f,TOp 
(where nip is given above and = ilb/f^m = 0.1429). 
We analyze clusters at the present-day epoch as well as 
their progenitors at higher redshifts. 

We repeated each cluster simulation with and without 
radiative cooling. The first set of "adiabatic" simulations 
have included only the standard gasdynamics for the 
baryonic component without dissipation and star forma- 
tion. The second set of simulations included gasdynamics 
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TABLE 1 

Simulated cluster sample at 



Name 


itvir 


Mvir 


A3l80m 


M200C 


iV/soOc 


^25000 


H 
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(^x ) 
















keV 


keV 


CLl 


1.911 


8.21 


9.62 


6.73 


5.31 


2.78 


6.9 


8.6 


CL2 


1.205 


2.06 


2.32 


1.75 


1.34 


0.75 


2.7 


2.6 


CL3 


1.227 


2.17 


2.45 


1.87 


1.32 


0.69 


2.7 


3.8 


CL4 


1.187 


1.97 


2.37 


1.65 


1.24 


0.63 


2.8 


3.2 


CL5 


1.016 


1.23 


1.41 


1.01 


0.78 


0.39 


1.8 


1.7 


CL6 


0.906 


0.87 


0.99 


0.78 


0.62 


0.34 


1.6 


2.5 


CL7 


0.947 


1.00 


1.22 


0.66 


0.48 


0.24 


1.4 


1.0 


CL8 


0.966 


1.06 


1.24 


0.88 


0.65 


0.34 


1.7 


2.1 


CL9 


0.795 


0.59 


0.68 


0.51 


0.35 


0.13 


0.8 


1.2 



and several physical processes critical to various aspects 
of galaxy formation: star formation, metal enrichment 
and thermal feedback due to the supernovae type II and 
type la, self-consistent advection of metals, metallicity 
dependent radiative cooling an d UV heating due to cos - 
mological ionizing background ijHaardt fc Madaulll996j) . 
We will use labels 'ad' and 'csf for the adiabatic simu- 
lations and simulations with cooling and star formation, 
respectively. The cooling and heating rates take into ac- 
count Compton heating and cooling of plasma, UV heat- 
ing, atomic and molecular cooling and are tabulated for 
the temperature range 10^ < T < 10^ K and a grid of 
metallicities, and UV in tensities using the Cloudy code 
(ver. 96b4. iFerland et al.. .199&) . The Cloudy cooling 
and heating rates take into account metallicity of the 
gas, which is calculated self-consistently in the simula- 
tion, so that the local cooling rates depend on the local 
metallicity of the gas. 

Star formation in these simulations was done using th e 
observationally-motivated recipe fe.g.. lKennicuttill998j) : 
P* — Pgas/**> 'with = 4 X 10^ yrs. Stars are allowed 
to form in regions with temperature T < 2 x 10^ K and 
gas density n > 0.1 cm~^. No other criteria (like the 
collapse condition V • v < 0) are used. We have com- 
pared runs where star formation was allowed to proceed 
in regions different from our fiducial runs. We consid- 
ered thresholds for star formation of n = 10, 1, 0.1, and 
0.01 cm~^. We find that thresholds affects gas fractions 
at small radii, r/r^ir < 0.1, but the differences are neg- 
ligible at the radii we consider in this study. The effect 
on the baryon fractions is small at all radii. 

Algorithmically, star formation events are assumed to 
occur once every global time step Ato ^ 10^ yrs, the 
value close to the o bserved star formation timescales 
fe.g.. lHartmannl2002() . Collisionless stellar particles with 
mass m, = Aio are formed in every unrefined mesh cell 
that satisfies criteria for star formation during star for- 
mation events. The mass of stellar particles is restricted 
to be larger than min(5 x 10^ h^^ M0, 2/3 x mgas), where 
"T-gas is gas mass in the star forming cell. This is done 
in order to keep the number of stellar particles compu- 
tationally tractable, while avoiding sudden dramatic de- 
crease of the local gas density. In the simulations ana- 
lyzed here, the masses of stellar particles formed by this 
algorithm range from « 10^ to 7 x lO^h^^ M©. 



Once formed, each stellar particle is treated as a single- 
age stellar population and its feedback on the surround- 
ing gas is implemented accordingly. The feedback here is 
meant in a broad sense and includes injection of energy 
and heavy elements (metals) via stellar winds and su- 
pernovae and secular mass loss. Specifically, in the sim- 
ulations analyzed here, we assumed that stellar initial 
mass function (IMF) is described by the Miller & Scalfll 
(jl979J functional form with stellar masses in the range 
0.1 - 100 Mq. All stars more massive than M, > 8 M© 
deposit 2 x 10^^ ergs of thermal energy in their par- 
ent cell^ and fraction /z = min(0.2, O.OlAf* - 0.06) of 
their mas s as metals, which crudel y approximates the 
results of iWooslev fc WeaveH l)199,^ . In addition, stel- 
lar particles return a fraction of their mass and met- 
als to the surrounding gas at a secular rate rhioss — 
Co(t-th\rth+Tn)- ^ with Co = 0.05 and Tq = 5 Myr 
(jjungwiert et al."2001). The code also accounts for the 
SNIa feedback assuming a rate that slowly increases with 
time and broadly peaks at the population age of 1 Gyr. 
We assume that a fraction of 1.5 x 10~^ of mass in stars 
between 3 and 8 M© explodes as SNIa over the entire 
population history and each SNIa dumps 2 x 10^""^ ergs 
of thermal energy and ejects 1.3 Mq of metals into the 
parent cell. For the assumed IMF, 75 SNII (instantly) 
and 11 SNIa (over several billion years) are produced by 
a lO'' M0 stellar particle. 

Throughout this paper we use estimates of gas and 
baryon fractions within different commonly used radii, 
defined by the total matter overdensity they enclose. We 
will use radii r25oo, '''500 ; '''200 enclosing overdensities of 
2500, 500, and 200 with respect to the critical density, 
Pcrit, as well as radii rigo and Vvh- enclosing overdensities 
of 180 and Avir with respect to the mean density of the 
universe. The latter is equal to Avir ~ 334 at z = 
and w 200 at z = 1 for the cosmology adopted in our 
simulations. 

The virial radius and masses of clusters within differ- 
ent radii are given in Table ^ For reference, we also 
give average emission- weighted temperature in the radial 
range between 80 kpc to rsoo of each cluster in adiabatic 
and cooling runs. The temperature is calculated as an 

^ No delay of cooling was introduced in these cells after SN 
energy release. 
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Fig. 1. — Comparison of the baryon fraction profiles in the adi- 
abatic simulations of the eight clusters CL2-9 done with the ART 
{solid lines) and entropy-conserving Gadget (long-dashed lines) 
codes. The gas fractions are normalized to the universal baryon 
fraction of each simulation, while radii of all clusters are normal- 
ized to their respective radius enclosing the overdensity of 200 with 
respect to the critical density. The bottom panel shows differen- 
tial gas fraction profiles, while the top panel shows corresponding 
cumulative profiles. The shaded bands represent la rms scatter 
around the mean for the eight clusters. The vertical arrows in the 
top panel show the radii enclosing overdensities of 2500, 500 (with 
respect to Pcrit): the virial overdensity and the overdensity of 180 
with respect to the mean density. 



emission-weighted average convolved with the Chandra 
energy response in the 0.5—7 keV energy band. 

3. COMPARISON OF THE ART AND GADGET 
SIMULATIONS 

Before we begin, it is important to discuss the de- 
gree of numerical uncertainty in the baryon fractions de- 
rived from cosmological simulations. We will compare 
simulations of eight lower-mass clusters (CL2-CL9) from 
the Table n performed with the eulerian AMR Adaptive 
Refinement T ree code and the entropy conserving SPH 
Gadget code iSoringel et al.ll200ll ISorineel fc Hernauisti 
[2002) performed in the adiabatic regime. The ART simu- 
lations are descr ibed above, while the G adget simulations 
are described bv lAscasibar et alJ l)200rif >. We compare the 
CL2-CL9 for which the SPH simulations are available. 
In each case the simulations have been started from the 
identical initial conditions and ran in the same cosmology 
with comparable spatial and mass resolutions, although 
Gadget simulations used somewhat different values of ilb 
for each individual cluster simulation. The number of 



dark matter particles is the same in the ART and Gad- 
get simulations, while the number of gas particles varies 
but is of order million for all clusters. Both gravita- 
tional and the minimum SPH softening are in the range 
e = 2 — 5h~^ kpc, depending on the number of particles 
within the virial radius. 

The average cumulative and differential gas profiles in 
the two sets of simulations are compared in Figure ^ 
The figure shows a rather remarkable agreement (to a 
couple percent) in the differential profiles at r > 0.4rvir, 
while at smaller radii there are systematic differences. 
The gas fraction profiles in the Gadget simulations are 
more centrally concentrated at r < 0.2rvir compared to 
the ART simulations. At the same time, the gas frac- 
tion is systematically higher in the ART simulations at 
r 0.2 — 0.4rvir. Comparison of the cumulative gas frac- 
tion profiles shows that despite the fact that the ART 
profiles are less centrally concentrated at small radii, the 
gas fraction within 0.4rvir is higher than in the Gadget 
simulations by « 5%. At larger radii, r > r25ooc, the 
systematic difference persists, although it is only about 
3 — 5%. Note also that the scatter about the mean is simi- 
lar in the ART and Gadget simulations. We checked that 
the differences in the mean profiles of the eight clusters 
also exist between profiles of individual clusters. They 
are therefore systematic and statistically significant. 

The difference is significantly smaller than the system- 
atic difference in gas fractions of ~ 10% between the 
eulerian and SPH codes, ob served in the Sant a Barbara 
cluster comparison project (|Frenk et al.lll999l) . Indeed, 
agreement in gas and dark matter density profiles at this 
level is quite remarkable. Nevertheless, it is important to 
keep these systematic differences in mind, both in gaug- 
ing implications of the results presented below and in 
observational analyses of gas fractions for cosmological 
constraints. As we argue below, this systematic differ- 
ence is smaller than the uncertainty in the simulations 
with cooling and star formation, associated with the un- 
certain fraction of baryons in stars. 

4. EFFECTS OF RESOLUTION 

Given that in these simulations we are trying to re- 
solve galaxy formation on kiloparsec scales while follow- 
ing formation of clusters self-consistently on scales of 
megaparsecs, it is reasonable to ask how the finite dy- 
namic range of the simulations affects our results. To 
this end, we have re-run one of the clusters used in our 
study (CL3) with higher resolution. The re-simulations 
use more agressive refinement criteria, which results in a 
factor of ~ 2.5 more refinement cells and a factor of two 
higher spatial resolution compared to the run used in this 
study. For example, the higher-resolution simulations at 
z = has « 2 X 10® grid cells on the 9th refinement level 
(comoving cell size of 1.22/i~^ kpc), while the runs used 
in this study did not have refinement cells beyond the 8th 
level. Gas fraction profile in the adiabatic re-simulation 
with such higher resolution agrees with the lower resolu- 
tion simulation to within ~ 1% at all resolved radii. We 
find a similarly small difference for the higher-resolution 
re-simulation of this cluster with the Gadget code with 
eight times more particles and softening of 0.5h~^ kpc 
(compared to the 2 — 5h^^ kpc in the simulations we use 
for comparison). In other words, it appears that for each 
code convergence in adiabatic runs is reached, so that dif- 
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ferent resolution does not explain systematic differences 
in the gas fraction profiles discussed above. 

The higher-resolution re-simulation of CL3 with cool- 
ing and star formation results in total baryon fraction 
profile that is within ~ 1 — 2% of that in the low- 
resolution simulation. However, we find that the gas frac- 
tion is ~ 3 — 5% higher in the higher resolution run. It is 
not clear what causes the difference, although we think 
that it is most likely be due to the increased heating by 
stellar feedback in the higher-resolution simulation. The 
difference is relatively small and is considerably smaller 
than the magnitude of the effects we discuss. Neverthe- 
less, the effects of resolution in the runs with star forma- 
tion will have to be studied further with higher-resolution 
runs. 

5. RESULTS 

We will start discussion of our results by considering 
an example of a single typical cluster from our sample 
(CL3). The cluster has virial mass similar to that of 
the Virgo cluster and has experienced a nearly equal- 
mass major merger at z « 0.6. We have analyzed three 
simulations of this cluster started from the same initial 
conditions but run including different physical processes. 
The first run was done in the "adiabatic" regime with- 
out radiative gas cooling and star formation. The second 
run included cooling and star formation. The third run 
is identical to the second at z > 2. At z < 2, however, 
gas dissipation and star formation have been artificially 
turned off. This run should be considered not as a re- 
alistic cluster model, but as a useful intermediate case 
between adiabatic run and run with full dissipation. For 
example, the fraction of baryons turned into stars within 
the virial radius is /* = 0.35 in the simulation with cool- 
ing and /, = 0.15 in the simulation with no cooling at 
2 < 2 (/, = for the adiabatic simulation). 

Figure 121 shows the cumulative and differential profiles 
of baryon and gas fraction for the three re-simulations 
of CL3. The figure shows several key qualitative ef- 
fects of dissipation. The most salient difference is in 
the gas fraction. The simulation with cooling has the 
largest stellar fraction and, correspondingly, the small- 
est /gas(< J'vir)- The difference in cumulative gas frac- 
tion is approximately constant at r > 0.15rvir. However, 
the lower panel of Figure |21 shows that the difference in 
differential profile of /gas increases monotonically with 
decreasing radius from less than 10% at r > 0.4rvir to 
~ 30 — 40% at smaller radii. This is because the cen- 
tral cluster galaxy has a significant effect on the overall 
stellar fraction and cumulative profile. 

Comparison of the total baryon fraction profiles for 
the three runs reveals additional differences. The 
baryon fraction in the adiabatic simulation is signifi- 
cantly smaller than the universal value at r < 0.2rvir, 
while it is larger than universal in the simulations with 
cooling. The cumulative baryon fraction in the adiabatic 
simulation reaches the universal value well beyond the 
virial radius. In the simulations with cooling the baryon 
fraction is close to the universal at the virial radius. Note 
that in all three runs the baryon fractions are within 10% 
of the universal value at r > 0.4rvir. The baryon fraction 
profiles are similar in all three simulations outside the 
virial radius. 

The main conclusion from this comparison, which also 
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Fig. 2. — Effects of cooling on the cluster baryon and gas fraction 
profiles at z = 0. The figure shows profiles for the three resimu- 
lations of the Virgo-size cluster CL3. Top 'panel: the cumulative 
baryon and gas fractions within a given radius in the simulations 
with cooling {solid, thick lines show /gas + f* while thin lines show 
/gas), cooling only at z > 2 (long- dashed lines), and no cooling 
(or adiabatic, dot-dashed line). The radius is in units of the virial 
radius, defined as the radius enclosing overdensity of Avir = 334 
with respect to the mean density of the universe. The fractions are 
in units of the universal baryon fraction, Qi,/^lm- Bottom panel: 
the same as the top panel, but for the differential profiles. 



applies to the comparison of the mean profiles presented 
below, is that the cooling and star formation change not 
only the amplitude but also the shape of the gas and 
baryon fraction profiles in the simulated clusters. 

This is to be expected because cooling and star for- 
mation lead to redistribution of both baryons and dark 
matter. The gas condensation results in a contrac- 
tion of the dark matter halo fe.g.. iZeldovich et alJll98(l 
iBlumenthal et alJll986|) and steepening of th e DM den- 
sity profile at r < O.lrvir ijGnedin et al.ll20 0j). More sig- 
nificantly, dissipation and resulting star formation con- 
vert a fraction of gas into stars with the efficiency which 
may depend on the cluster radius. Dynamically, the stel- 
lar component is coUisionless and its evolution will differ 
from that of the gas. For instance, cluster galaxies may 
loose their gas via ram pressure stripping. The lost gas 
can then mix with the intracluster gas at large radii. 
The stellar components of galaxies, on the other hand, 
can sink to the center via dynamical friction and build 
up the massive central cluster galaxy. Likewise, the de- 
tailed dynamics and relaxation process of gas and stars 
during mergers is rather different. For example, the gas 
experiences heating via shocks, while stars participate in 
coUisionless violent relaxation. 
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Fig. 3. — Differential {bottom panels) and cumulative {top panels) baryon and gas fraction profiles for the nine clusters used in our 
analysis at z = {left column) and z = 1 {right column). The dot-dashed lines show the mean profiles in the adiabatic simulations averaged 
over eight clusters (CL2-CL9). The solid lines show the total mean baryon fraction (gas+stars) profile, while the long-dashed lines show 
the gas fraction profile in the simulations with cooling and star formation. The shaded bands show the Icr rms scatter around the mean 
for the eight clusters. The thin lines in the top panels show the corresponding profiles in the Coma-size cluster simulation (CLl). These 
profiles have not been used in the averages due to systematically different gas and baryon fractions. The vertical arrows in the top panels 
show the radii enclosing overdensities of 2500, 500, 200 (with respect to Pcrit)i and the overdensity of 180 with respect to the mean density. 



Figure|21shows the comparison of the mean baryon and 
gas fraction profiles at z = 1 and z — for our entire 
sample of clusters simulated in the adiabatic regime and 
with dissipation. The profile for the Coma-size cluster 
(CLl) is plotted separately because in this cluster stel- 
lar fraction (and hence the gas fraction) is significantly 
different from the rest of the smaller-mass clusters, as 
we discuss below. However, the difference is mainly in 
the amplitude of the profile. Qualitatively, the profiles 
of CLl are similar to the mean profile for the other eight 
clusters. 

The trends in the mean profiles are similar to those 
seen in Figure |21 The cumulative baryon fraction profiles 
are fiat at r > 0.4rvir. At these radii /bar is within « 5% 
of the universal value for both adiabatic and CSF runs, 
although in the former the /bar is systematically below 
the universal value, while in the latter it is systematically 



above. In the cores of the low-mass clusters, the cumu- 
lative baryon fraction is considerably below (above) the 
universal value in the adiabatic (CSF) runs. This shows 
that dissipation and star formation result in a substan- 
tial redistribution of baryons within the virial radius. In 
the Coma-size cluster the baryon fraction is still close 
to the universal value within ^2500, but at smaller radii 
the trends are similar. Note that the differential baryon 
fraction in the CSF simulations is systematically below 
the universal value by w 10 — 20%. The profile, how- 
ever, is rather flat at r > ^2500- Both the cumulative and 
differential gas fractions increase monotonically with ra- 
dius. If a similar trend exists in real clusters this implies 
that the correction of gas fraction for stars in general de- 
pends on the cluster-centric radius and a single assumed 
value of /, for all radii may not give the universal baryon 
fraction. 
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Fig. 4. — The total baryon and gas fractions of individual clusters in units of the universal value within different radii (clockwise from 
the top left panel: r < r2500i ^ < ^"500, ^ < ^vir, in the annulus r2500 — ^50o) as a function of the cluster virial mass. Solid triangles show 
the gas fractions in the adiabatic simulations. The solid circles show the gas fraction, while open circles show the total baryon fraction in 
the simulations with cooling and star formation. Solid and open squares show the gas and baryon fractions in the resimulation of CL3, in 
which cooling was turned off at 2 < 2. Note the systematic trend with mass for the gas fractions in the simulation with cooling. 



Figure|31also shows that the variance around the mean 
for the cumulative profiles is rather small, especially at 
large radii. The scatter is somewhat larger at z = 1, 
probably because the clusters are on average less relaxed 
at the earlier time. However, there is no significant qual- 
itative difference between the profiles at two different 
epochs. 

Figure 0] shows the baryon and gas fractions as a func- 
tion of the cluster virial mass within r25oo , ''500 j E^nd Tvir 
and within the radial range [r25oo — ''500] • The values of 
the baryon and gas fractions plotted in this figure are 
listed in Table El The fractions in the adiabatic simula- 
tions are marked 'ad', while those in the simulations with 
cooling and star formation are marked 'csf (e.g., /g^s). 
In the table we also provide the mean values and scat- 
ter for our cluster sample. Note, however, that these are 



averages over clusters of different masses and should be 
interpreted with caution given the trends of the baryon 
and gas fractions with cluster mass discussed below. 

In adiabatic simulations, the baryon fraction appears 
to be independent of mass for all radial ranges. It is 
systematically below the universal value by « 15% and 
~ 5% within r25oo and rsoo, respectively. For r < r^ir, 
the baryon fraction is within 5% of the universal value. 
Note also that in the adiabatic runs, the baryon fraction 
within the radial range of r25oo — '"500 is close to the 
universal value on average, even though the cumulative 
fractions within these radii are systematically below. 

In the simulations with cooling and star formation the 
baryon and gas fractions show (weak) trend with mass. 
The gas fraction appears to systematically increase with 
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Fig. 5. — Evolution of the baryon and gas fractions within different radii for the eight clusters CL2-CL9 from z = 1 to z = 0. The radii 
used are the same as in Fig.|^ The long-dashed lines shows the gas fraction in the adiabatic simulations. The dot-dashed and solid lines 
show the gas and baryon fraction profiles in the simulation with cooling and star formation. The gas and baryon fractions are normalized 
to the universal fraction of baryons assumed in the simulations. The figure shows that the evolution of baryon fraction at these redshifts 
is weak. Note that each cluster at least doubles its mass over this redshift interval. 



increasing cluster mass as 

/gas = 0.21og(Mvi,/10i4/i-i Mo) + /o, (1) 

where fo is « 0.4, 0.5, 0.6, 0.65 for r < r25oo, ''soo, r^h , 
and r25oo ~ ^500? respectively. The significance of this 
trend and the value of the slope are uncertain given the 
small size of the sample. The baryon fraction in the 
CSF runs is larger than the universal value within r^oo, 
but it is within ?» 2 — 3% of the universal value within 
Tvir- The baryon fraction, /bar(< ''vir), in these runs is 
also not biased low like the corresponding fraction in the 
adiabatic runs but gives the average very close to the 
universal. 

The differential fractions within the radial range 
'"2500 ~ ''soo behave somewhat differently because they 
are not affected by the concentration of baryons in the 
cluster center. In this radial range, the baryon frac- 
tion in the adiabatic runs is approximately unbiased. 



In the CSF runs the baryon fraction is systematically 
below the universal value but its mass dependence for 
M > lO^^ft.^^ M0 seems to be weaker than for the cu- 
mulative fraction within similar radii. In all panels open 
and solid squares show the baryon and gas fractions for 
the re-simulation of CL3 with cooling turned off at z < 2. 
As could be expected, the fractions in this case are inter- 
mediate between the adiabatic simulation and run with 
cooling and star formation. 

Finally, in Figure [S] we show evolution of the average 
baryon and gas fractions from z = 1 to z = within 
radii shown in the previous figure in the adiabatic and 
CSF runs of clusters CL2-CL9. The figure shows that 
the total baryon fractions are quite stable in this redshift 
interval. This is remarkable given that most clusters un- 
dergo at least one major merger at z < 1. There is a 
weak decrease of the gas fraction with time in the CSF 
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simulations, as some of the gas is cooling and is turned 
into stars. The evolution of /gas is « 10 — 20% from z = 1 
to z = 0. 

6. COMPARISON TO PREVIOUS WORK 

Given the systematic differences in gas fractions in the 
Santa Barbara cluster comparison l|Frenk et al.lll999H . it 
is interesting to compare the results presented in the 
previous section to other recent studies to assess the 
current state of affairs. The comparison is relatively 
straightforward for the adiabatic simulations because all 
codes should simulate evolution using the same equa- 
tions of gas and coUisionless dynamics. Any differences 
can therefore be attribu ted to the differe nce in the ac- 
tual numerical scheme. lEke et aTl (|,1998) studied SPH 
simulations of ten massive clusters formed in the concor- 
dance ACDM cosmology. They found gas fractions of 
« 0.99, 0.87, and 0.83 within 3rvir, ?'vir, and O.Srvir (note 
that r5oo ~ 0.5rvir), respectively, which evolve only very 
weakly from z = 1 to z = 0. These numbers are simi- 
lar to those o btained w ith a number o f other SPH codes 
llFrenk et all [1999: Bia lek et all 120011: iMuanwong et all 
We also find no detectable evolution of baryon 
fraction in adiabatic simulations over this period of time. 
However, the gas fractions in our simulations are larger: 
« 1.00 at r > Brvir, and 0.97 ± 0.03 and 0.94 ± 0.03 
for r < Tvir and r < r^oQ. Gas fractions in simulations 
with the entropy-conserving Gadget co de (see Fig . JTjl ar e 
larger by about 5% than those found bv lEke et al.ll|1998|) : 
~ 0.92 at r < Tvir and « 0.88 if measured within O.Srvir. 
This is consistent with the adiabatic simulations (run 
with the entropy conserving version of the G adget) pre- 
sented in a recent study bv lKav et al.l l)2004f) . who find 
the average gas fraction of ~ 90% of the universal value 
within the virial radius. The average gas fraction profile 
for the adiabatic simulations in their Figure 6 is in good 
agreement with the Gadget simulations presented here. 

The studies of baryon fraction in simulations with cool- 
ing and star formation are fewer and the source of the 
differences is usually not as clear, because it can be 
attributed both to a different implementation of these 
processes and to th e difference betwe e n num erical gas- 
dynamics schemes. IMuanwong et al.l l)2002j) compared 
the baryon and gas fractions as a function of cluster 
virial mass in adiabatic simulations and re-simulations 
in which gas was allowed to cool or was pre-heated (see 
their Fig. 3). For the mass r ange of our cluster sam- 
ple, in simulations with coohng IMuanwong et alJ 1)20021 ) 
obtain baryon fractions of « 0.85 — 0.9 of the universal 
value within the virial radius. This is slightly larger than 
in their adiabatic simulations, but is considerably lower 
than in our simulations with cooling and star formation 
(1.02±0.02, see Tahle^. iValdarninil l )2003|) compares gas 
fractions in in TREESPH simulations that include cool- 
ing and star formation to the observational estimates of 
Ernaud & Evrard (1999) and finds a reasonable agree- 
ment. iValdarninil (|2003J also finds indications of a trend 
of increasing /gas with increasing cluster mass similar to 
the trend we find in our simulations. The actual gas frac- 
tions in this study are, however, somewhat lower than 
our values, which may be due to a lower assumed univer- 
sal baryon fractions {il\^h^ = 0.01 5 — 0.019) c ompared to 
the value assumed in this study. iKav et ahl |2004) ana- 
lyze a sample of 15 clusters simulated using the entropy- 



conserving Gadget code with cooling, star formation, and 
feedback. Their average gas fraction profile for these 
simulations agree well with the gas fraction profile of 
our small- mass clusters (CL2-CL9), but is considerably 
lower th an that for our Co ma-size cluster (CLl). Most 
recently, lEttori et al.l l|2004D presented a similar analysis 
for the simulations with cooling and star formation of a 
statistical sample of clusters simulated with the entropy- 
conserving version of the Gadget code. These authors, 
quote total baryon fractions of /bai (< ''vir) ~ 0.93 — 0.95 
of the universal value in their simulations. This is closer 
to the values in our simulations. Their gas fractions are 
0.75 — 0.8 — higher than in our clusters due to a different 
implementation of cooling and star formation. 

These comparisons and the comparison with the Gad- 
get code presented in §|3indicate that there are still sys- 
tematic differences between the eulerian and SPH codes, 
with the latter predicting systematically smaller baryon 
fractions. The difference appears to be the smallest 
when eulerian simulations are compared to the entropy- 
conserving Gadget SPH code. 

7. DISCUSSION AND CONCLUSIONS 

We have presented results of the analysis of baryon 
and gas fractions in the high-resolution simulations of 
nine galaxy clusters formed in the concordance AGDM 
cosmology. The clusters span the range of virial masses 
from Mvir = 6 X 10" to 8 x lO^/i'^ Mq and have not 
been selected to be in the highly relaxed state at the 
present epoch. We study the effects of cooling and star 
formation by comparing simulations of individual clus- 
ters done with and without these processes included. 

We show that radiative cooling of gas and subse- 
quent star formation significantly modify distribution of 
baryons within the virial radii of galaxy clusters in our 
simulations. The effect is twofold. First, cooling results 
in a cond ensation of gas and da r k matter in the centers 
of halos JZddlTvich et. a].l[T98?i iBarnes fc Whitfl IT981 
Bluinenthal et alJll986HJesseit et al.ll2002HGnedin et alJ 
20041) . As the gas cools in the cluster progenitor, the 
baryons condense near the cluster center (the central 
cluster galaxy). The condensing gas has to be replaced 
by the gas from larger radii. The net result is that the 
baryon fraction is affected at radii up to ^ rvh-. For 
example. Figure |21 shows that the cumulative baryon 
fraction in the simulations with cooling is systematically 
higher than in the adiabatic runs at r < 1.5rvir with 
the difference of « 5% at rvir. Dissipation thus allows 
clusters to collect a larger fraction of baryons within 
their virialized region. Remarkably, in these simulations 
baryon fraction within the virial radius is very close (to 
w 2 — 3%) to the universal value. Note that for the dif- 
ferential baryon fraction, fhniir), the effect at large radii 
is opposite. In the simulations with coohng /bar('') is 
systematically smaller than in the adiabatic simulations 
at r > O.lrvir- The overall increase in the cumulative 
baryon fraction is therefore driven by the baryons in the 
central cluster galaxy. 

Second, a fraction of gas is converted into coUisionless 
stellar component. The ratio of stellar to gas mass de- 
pends on the distance to the cluster center. In cluster 
cores (r < 0.2rvir) stellar fraction dominates the baryon 
budget, /* > 50 — 70% depending on cluster mass. At 
larger radii most baryons remain in the form of gas with 
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TABLE 2 

Baryon and gas fractions in cluster simulations at the present day epoch 



Cluster 


r < r2500 

rad pcsi 
J bar J gas 


csf 
.^bar 


ad 
^ bar 


>" < '"500 
/gas 


csf 
/ bar 


^2500 - ^5 
rad pcs{ 
Jbar Jgas 


00 ^ 

/ bar 


ad 
/ bar 


r < rvir 

/gas 


est 
/bar 


CLl 


0.86 


0.56 


1.06 


0.98 


0.69 


1.02 


1.10 


0.87 


0.98 


1.04 


0.76 


1.05 


CL2 


0.73 


0.44 


1.17 


0.88 


0.57 


1.05 


1.08 


0.67 


0.79 


0.96 


0.63 


1.01 


CL3 


0.95 


0.49 


1.18 


0.96 


0.57 


1.05 


0.97 


0.71 


0.83 


0.96 


0.64 


1.01 


CL4 


0.86 


0.50 


1.18 


0.92 


0.60 


1.02 


0.99 


0.74 


0.80 


0.95 


0.68 


0.99 


CL5 


0.90 


0.48 


1.16 


0.95 


0.62 


1.06 


0.99 


0.72 


0.99 


0.97 


0.68 


1.02 


CL6 


0.92 


0.34 


1.26 


0.95 


0.45 


1.06 


0.98 


0.63 


0.73 


0.97 


0.55 


1.02 


CL7 


0.87 


0.45 


1.18 


0.92 


0.60 


1.08 


0.98 


0.82 


0.92 


0.99 


0.68 


1.04 


CL8 


0.83 


0.47 


1.29 


0.97 


0.56 


1.07 


1.12 


0.68 


0.82 


0.97 


0.64 


1.02 


CL9 


0.72 


0.41 


1.45 


0.91 


0.53 


1.14 


1.03 


0.67 


0.79 


0.94 


0.61 


1.00 


mean 


0.85 


0.46 


1.22 


0.94 


0.58 


1.06 


1.02 


0.72 


0.87 


0.97 


0.65 


1.02 


scatter 


0.08 


0.06 


0.11 


0.03 


0.07 


0.03 


0.06 


0.08 


0.09 


0.03 


0.06 


0.02 



/* < 10 — 20% at r > r25oo- The collisionless dynamics 
of stars is different from the dynamics of gas in detail. 
For example, the gas in a galaxy or a group can be ram- 
pressure stripped and deposited at larger radii, while the 
stars can sink to the center via dynamical friction. Dur- 
ing mergers, the gas is shock-heated, while stellar par- 
ticles relax via the violent relaxation. When stars are 
formed, a smaller fraction of baryons is subject to strong 
shocks that may prevent accretion of some of the gas into 
the inner regions of clusters. The differences in dynamics 
can modify the total baryon fraction profile and the ratio 
of stellar to gas mass at different radii. 

To gauge how robust are the predictions of numerical 
simulations, we have compared the gas fraction profiles 
for the eight clusters in our sample, simulated in the adi- 
abatic regime with the entropy-conserving Gadget and 
the ART codes. The profiles in the two sets of simula- 
tions agree to better than « 3% outside the cluster core 
{r/rvir ^ 0.2), but differ by up to 10% at small radii. 
The main difference is thus in the cluster core where the 
gas profiles in the Gadget simulations are more centrally 
concentrated and the gas fractions are correspondingly 
larger than in the ART runs. This is reflected in system- 
atic difference in the cumulative baryon fraction profiles 
at r > ^2500: 95% in the ART runs compared to 92% in 
the Gadget simulations. The discrepancy is smaller than 
the difference in gas fractions of ~ 10% between SPH and 
eulerian codes in t he Santa Barbara comparison project 
l)Frenk et al."1999^. Nevertheless, it is systematic (simi- 
lar for all eight compared clusters) and one has to keep it 
in mind when interpreting or using results of numerical 
simulations. Note also that this comparison was done 
for the simplest case of the adiabatic gas dynamics. The 
differences in baryon and gas fractions are likely to be 
larger when cooling and star formation are included. We 
discussed how our results compare with some of the re- 
cent simulation results in the previous section. We also 
plan to compare our results with the results of the Gad- 
get simulations with cooling and star formation in the 
near future. 

One should note that cluster simulations may suffer 
from the "overcooling problem." Indeed, the fraction of 



baryons in the cold gas and stars within the virial ra- 
dius of clusters at z = in our simulations is in the 
range ^ 0.25 — 0.35, at least a factor of two higher than 
observational meas ureme nts for the systems of the mass 
range we consider ()Lin et al, 2004) . Such high fractions 
of baryons in condensed cold gas or stars are generic re- 
sults of cosmological simulations, although specific nu- 
merical val ues vary (e.g ., Suainohara Ik Ostrikcr 199S 
Lewis et al .' 2000; Pear ce et all I2000t iDave et all 12003 
Ettori et al . 2004). We have done an extensive conver- 



gence study, varying the numerical resolution, implemen- 
tation of cooling, star formation and stellar feedback, and 
assumed baryon fraction. However, we have not been 
able to reduce stellar fractions considerably below the 
values prese nted in this s tudy. 

Note that iMotl et all l)2004ll have recently presented 
eulerian AMR simulations of two massive galaxy clus- 
ters (Mvir « 2 X 10^^/1"^ Mq) with the fraction of cold 
(T < 15000 K) condensed gas of « 10%. It is not yet 
clear whether this represents a discrepancy with our re- 
sults. The number of simulated massive clusters is clearly 
small. At the same time, the two simulated clusters an- 
alyzed by IMotl et a l. ( 2004) have mass a factor of two 
larger than our most massive cluster. If the trend of in- 
creasing /gas (and correspondingly decreasing f.^,) with 
increasing cluster mass (see Fig. 0J) continues to larger 
masses, equation gives /p^a.s y 90% for the cluster 
masses considered bv IMotl et al.l l)2004D . in agreement 
with their results. 

The efficient dissipation and persistence of late star for- 
mation may indicate that some mechanism suppressing 
gas cooling is needed. The problem appears to be in the 
central brightest cluster galaxy, which contains a larger 
fraction of cluster mass than is measured in observations. 
The rest of the cluster galaxies have stellar masses in rea- 
sonable agreement with observations (Nagai & Kravtsov, 
in preparation) . The overcooling problem in simulations 
may thus be related to the puzzli ng absence of cold ga s 
in the observed cluster cores ('e.g..lPeter son et al.ll2003D . 
At the same time, some fraction of cluster stellar mass 
in low-surface brightness intracluster component may be 
missed in observations, which could reduce the apparent 
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discr epancy of observed f^ , with the sim ulation resuhs 
Ce.g.. IConzalez et a,1.ll2nnl IT.in MohHl2nn4ll . For ex- 
ample, baryon fractions within virial radii of rich clus- 
ters do not seem to add up to the unive rsal value from 
the WMAP observations (Spcrgcl ct al.'2003), as they do 
in both adia batic and dissipative simulations (see, e.g., 
IEttorill2f)0l . 

Interestingly, gas fractions measured in our CSF sim- 
ulations are consist ent with the observed ga s frac- 
tions in clus ters (Evrard 1997; Mohr et al. 1999"; 
lArnaud fc Evra rd 1999: Roussel ct al. 2000). For exam- 
ple. lEvrardl ljl997') presents estimates of the gas frac- 
tions within radius of rsoo for a sample of galaxy clus- 
ters spanning a wide range of the ICM temperatures. 
The gas fractions vary from ~ 50 to 85% of the univer- 
sal value (assuming fib = 0-043 and h = 0.7) in rough 
agreement with our results. The values of gas fraction in 
units of the universal value, /bar = 0.143, est i mated from 
the virial scaling relations by IMohr et al.l ljl999(l are^ 
/gas(< ''50o)//bar ~ 0.5 — 0.7 for clustcrs with tempera- 
tures of Tx « 4 keV and /gas(< r5oo)//bar - 0.65 - 0.85 
for Tx ~ 8 keV. The latter is also consistent with the 
measurements based on the SZ cluster obse rvations and 
X-ray temperature of Tx 5 keV clus t ers llGrego et al.l 
I200I . More recently, ISanderson et al.l 1)20031) presented 
analysis of gas fractions in a large sample of clusters span- 
ning a wide range of masses. The resulting values of gas 
fraction depend on the assumptions made, but overall the 
values are larger than in our simulations. For instance 
the mean gas fraction for their clusters is 0.13 ± 0.01 
or in units of our universal baryon fraction 91%. Simi- 
larly to our simulated clusters, the real clusters in their 
sample show gas fract ions monotonical l y incr easing with 
cluster-centric radius. ISanderson et alJ i|2003|) also found 
a trend of increasing gas fraction with increasing cluster 
mass (see also Mohr et al. 1999; Arnaud & Evrard 1999) 
roughly consistent with the trend for our sample of simu- 
lated clusters (see Figure 0] and eq. pOj). A caveat is that 
gas fractions in these studies have been estimated either 
using the /3-model fits to the data or virial scaling rela- 
tions rather than directly measuring gas and total mass 
within rsoQ. 

We also note that IBrvanI l)2000|) used observational 
data compiled from the literature to argue that there 
is a trend of decreasing stellar mass fractions with in- 
creasing cluster mass: /, = 0.042(r/10 keV)~°-^^, where 
T is the emission-weighted ICM temperature. Given 
that T oc M^/'^ according to the expected virial scal- 
ing, the corresponding scaling as a function of cluster 
mass, /* oc Af~°-^'^, and its normalization are in good 
agreement with our simulations. 

More recently, high sensitivity of the Chandra and 
XMM-Newton satellites allowed reliable measurements 
of temperature profiles and direct estimates of gas frac- 
tions from the full hydrostatic equilibrium analyses. For 
instance, Allen et al. (2004) present gas fraction pro- 
files at r < r25oo for a sample of 26 massive clusters 
{Tx > 5 keV) at different redshifts (0 < z < 1). They 
quote gas fractions /gas(< ?'250o) ~ 0.12 for the cosmol- 
ogy adopted in our study. The corresponding gas frac- 
tion for the massive cluster (CLl) in our study is 0.08. 
This lower value may be due to the fact that most of 

We have calculated the gas fractions for h = 0.7. 



the clusters in the lAllen et al.l l)2004j) are more massive 
than CLl^ or may indicate overcooling in simulations. 
It may also be simply due to the higher value of the 
universal baryon fraction in the real Universe. For ex- 
ample, if we us e the value preferred by the WMAP data 
(|Spergel et al.ll2003,) . n^/il^ ~ 0.17, the CLl gas frac- 
tion in units of the universal value agrees with the ob- 
served values in lAllen et al.l l)2004|) . A larger sample of 
massive simulated clusters or comparison with a sample 
of lower mass clusters is needed for a reliable conclusion. 

The results presented in this paper have a number of 
implications for the analyses that use clusters to obtain 
cosmological constraints. First of all, our results show 
that cooling changes the cumulative baryon and gas frac- 
tions appreciably even at the virial radius. Within the in- 
ner regions of clusters (r < »'25oo) the difference between 
cooling and adiabatic simulations of the same cluster can 
be as large as ~ 20 — 40%. The effects are smaller for 
more massive clusters or if the real clusters have lower 
fraction of condensed baryons than clusters in our simu- 
lations. However, even for the simulation in which stellar 
fraction is twice smaller (/* w 0.15 in the re-simulation of 
CL3 with cooling turned off at z < 2) we find 20% change 
in gas fraction and « 7% change in the total baryon frac- 
tion within r25oo compared to the adiabatic run. These 
differences are still large if baryon fractions are to be used 
for precision cosmological constraints. Our analysis also 
shows that in simulations with cooling the total baryon 
fraction can either be larger or smaller than the univer- 
sal fraction depending on the chosen radius and cluster 
mass. Therefore, there is no single universal correction 
for converting the measured baryon fraction to the uni- 
versal value. These results suggest that care should be 
taken in observational analyses to estimate baryon frac- 
tions, and the systematic errors are likely mitigated if the 
measurements are made at radii least affected by cooling 
and least sensitive to the uncertain contribution of stars 
to the baryon budget. 

Cooling and star formation decrease the fraction of 
hot intracluster gas within the virial radius of the clus- 
ter. This decrease may affect the Sunyaev-Zeldovich (SZ) 
fluxes and their relation to cluster mass. Upcoming large 
SZ surveys will have a large fraction of objects with 
masses of a few x 10^'' Mq . In such clusters star formation 
can decrease the gas fraction by ~ 20 — 40% depending 
on cluster mass. The trend of gas fraction with mass 
will also modify the slope of the relation between SZ flux 
and cluster mass. The trend is weak but may need to 
be taken into account to ensure accurate cosmological 
constraints. 

Given the importance of these issues for the use of clus- 
ters as cosmological probes, further effort from both the- 
orists and observers is needed to evaluate effects of cool- 
ing and to make robust measurements of baryon fractions 
(gas-|-stars) in clusters. Theoretically, it is important 
to simulate larger samples of clusters over a wide range 
of masses and with high resolution and to investigate 
the role of various processes (such as, for example, ther- 
mal c onduction and AGN feedback; e.g.. iBruggen et all 
I2OO5I) in suppressing cooling and star formation in clus- 
ter cores, as well as processes such as helium segrega- 

^ There are some clusters in their sample that have /gas(< 
^•2500) < 0.1 
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tion which can affect clu ster baryon fractions l)Qin fc Wul 
l2f)nntlChii7,hov T>oehl i2nn41. However, it IS even more 
important to first understand how large is the difference 
between the stellar fractions in observed and simulated 
clusters. If the fraction of the diffuse intr acluster stars 
is as large as 50% ( Gonzalez et al. 2003: iLin fc Moh^ 
12004 see however IZibetti et all l2005j) . L m our simu- 
lations actually agrees with observations. Careful obser- 
vational calibration of stellar masses in clusters (see, e.g., 
IVoevodkin et a l. ••2002: ' Lin et al.l 20041 using deep images 
in red or infrared bands would go a long way in clarifying 
this issue. 
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